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ABSTRACT 

We report on the serendipitous discovery of a 442-Hz pulsar during a Rossi 
X-ray Timing Explorer (RXTE) observation of the globular cluster NGC 6440. 
The oscillation is detected following a burst-like event which was decaying at the 
beginning of the observation. The time scale of the decay suggests we may have 
seen the tail-end of a long-duration burst. Low- mass X-ray binaries (LMXBs) 
are known to emit thermonuclear X-ray bursts that are sometimes modulated by 
the spin frequency of the star, the so called burst oscillations. The pulsations 
reported here are peculiar if interpreted as canonical burst oscillations. In par- 
ticular, the pulse train lasted for ~500 s, much longer than in standard burst 
oscillations. The signal was highly coherent and drifted down by ~2xlO~ 3 Hz, 
much smaller than the ~Hz drifts typically observed during normal bursts. The 
pulsations are reminiscent of those observed during the much more energetic "su- 
perbursts" , however, the temporal profile and the energetics of the burst suggest 
that it was not the tail end nor the precursor feature of a superb urst. It is possible 
that we caught the tail end of an outburst from a new 'intermittent" accreting 
X-ray millisecond puls ar, a phenomenon wh ich until now h as only been s een in 



HETE J1900. 1-2455 (jGallowav et all 120071 ). We note that Kaaret et all (120031 ) 



reported the discovery of a 409.7 Hz burst oscillation from SAX J1748.9— 2021, 
also located in NGC 6440. However, Chandra X-ray Observatory imaging indi- 
cates it contains several point-like X-ray sources, thus the 442 Hz object is likely 
a different source. 

Subject headings: — stars: neutron — X-rays: bursts — globular clusters: indi- 
vidual (NGC 6440) 
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Introduction 



The discovery of millisecond spin periods of neutron stars in Low Mass X-ray Binaries 
(LMXBs) with the Rossi X-ray Timing Explorer (RXTE) has helped elucidate the nature 
of these sources. Neutron star LMXBs consist of a neutron star accreting from a low mass 
companion. As material (mostly H and He ) is accreted onto the star and gets compressed, it 
eventually ignites and burns unstably (see IStrohmayer fc Bildstenl 120071 ). This phenomenon 
is observed a s a Type I X-ray Burst. Type I X-ray bursts have been observed from over ~70 



LMXBs (see iLiu et al.l 120061 . and references therein). The recurrence ti me of these bursts 



varie s but in some cases it can be as frequent as every few hours (see iGalloway fc Muno 
20071 . for examples). Occasionally it is possible to observe the spin of the neutron sta r 



modulating the burst emission - the so called "burst oscillatio ns" (IStrohmayer et al.lll996f) . 
Burst oscillations have been observed from ~18 LMXBs (see IGalloway fc Mund 120071 . for 
examples) . 

X-ray bursts from LMXBs have been discovered which are ~1000 times longer, and 
thus, that much more energetic, than canonical Type I X-ray bursts. They are aptly 
named "superbursts" . Supe r bursts are believed to occur by the unstable burning of c arbon 



(IStrohmayer fc Brown! 120021 ; ICumming fc Bildstenl 120011 ) . IStrohmayer fc Markwardtl (120021 ) 
discovered highly coherent pulsations during a superburst from 4U 1636—53. The pulse train 
lasted for ~800 s, as opposed to the ~10 s long pulse trains observed in Type I X-ray bursts. 
Superbursts are m uch rarer than Type I X-ray bursts. Thus far ten have been observed from 
eight LMXBs (see lin't Zand et all 12004 ). 



2. Observations 

Between 2005 March 7 and 2005 July 21 the PCA Galactic Bulge Scan Survejfl dis- 
covered an outburst from the direction of the Globular Cluster NGC 6440. A follow up 
1.8-ks long pointed RXTE observation (observation identification number 91050-03-07-00) 
was initiated on 2005 June 14. The data presented here were acquired from the Proportional 
Counter Array (PCA) on board RXTE. The PCA consists of five identical and independent 
proportional counter units (PCUs). Each PCU is a 90% Xenon/10% Methane gas filled 
proportional counter. Each PCU has a collimated l°xl° field-of-view, 256 spectral chan- 
nels in the 2-60 keV range, and a limiting temporal resolution of ~1 /is. Only two PCUs 
(PCU and PCU 2) were operational throughout the span of the observation. The data were 



4 http : / /lheawww . gsf c . nasa . gov/users/ craigm/galscan/ 
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taken in E_125us_64M_0_ls mode, which returns events to a limiting resolution of 125 /is and 
with moderate (64 channels as opposed to the full 256 channels) spectral resolution. This 
mode was used because it is not as susceptible to buffer overflows during high countrate 
data as compared to the less restrictive modes (e.g. GoodXenon). Using the pointing of the 
observation (17 h 48 m 52;8, -20°21'32") and the planetary ephemeris DE200, the times of raw 
events were corrected to the solar system barycenter. Binning the events into a 1-s resolution 
lightcurve reveals a quickly decaying burst (Fig. CD). 



3. Analysis and Results 



This burst was reminiscent of canonical Type I X-ray bursts from LMXBs, therefore we 
searched this event for burst oscillations. We rebinned the raw events using the full spectral 
band pass into a time series with 0.5/1024 s res olution, which yiel ds an equivalent Nyquist 
frequency of 1024 Hz. A Leahy normalized (see lLeahy et al.lll983l ) power density spectrum 
(PDS) is displayed in Figured] (inset). A prominent peak is clearly seen at 442 Hz. The 
probability of this peak occurring by chance after accounting for the number of trials (the 
total number of frequency bins in our PDS) is ~2 x 10 -9 . 



3.1. The NGC 6440 Field 



NGC 6440 harbors t he bright X-ray tra nsient SAX J1748. 9-2021, which exhibited a 
409.7 Hz burst oscillation (iKaaret et al.ll2003l ). However, there are many other X-ray sources 
in the cluster (see Fig[2]) . Chandra X-ray Observatory ( CXO) imaging by lPooley et al.l (120021 ) 
revealed 24 X-ray sources in NGC 6440, and they concluded that 4-5 of these sources are 
likely quiescent LMXBs. Thus, the phenomenon we discovered is most likely from a different 
object in the cluster. 



3.2. Timing Analysis 

To study the frequency evolution of the pulse we calculated a dynamic Z 2 statistic .The 
Z 2 statistic is analogous to the Fourier power spectrum, with the advantage that the data 
need not be binned, thus allowing us to oversample our data. The Z 2 statistic is defined as: 

9 Miarm ^7 

zi = — - y y | e i27rfa ^i 2 (i) 

1 k=Q i=0 
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Fig. 1. — Burst lightcurve binned with 1-s re solution and using the full PCA bandpass. 
Inset: "Leahy" normalized (ILeahy et al.l Il983l ) power spectrum of the entire observation. 
Notice the significant peak at 442 Hz. 
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Fig. 2. — Archival CXO image of NGC 6440. Notice that the cluster contains many X-ray 
point sources. 
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where iV 7 is the total number of photons in each interval, iVharm is the total number of 
harmonics that one deems significant, Uj is the frequency searched over, and £j is the event 
time. The factor in front of the summations normalizes the Z 2 statistics in an analogous 
way to Leahy normalizing a PDS. We calculated the dynamic Z 2 power spectrum using a 
200-s long window, which was translated across the observation with a step size of 16 s. Our 
dynamic Z 2 power spectrum is displayed in Fig. [31 Notice that the pulsations where highly 
significant for 576 s. As is common for canonical X-ray bursts, the frequency drifts; however, 
unlike the large ~l-2 Hz drifts seen in those, the pulsation here drifts down in frequency 
only by ~2.1xl0~ 3 Hz in 576 s. 

To further quantify this frequency evolution we determined how the phase of the pulsa- 
tions varies as a function of time. For each interval with Z\ > 16 we generated a pulse profile 
by folding the events in that interval on the frequency determined from the Z 2 statistic. We 
then cross correlated these pulse profiles with a sinusoid of fixed phase to determine the 
pulse phase as a function of time. We can model the pulse phase ((f)) at a given time (t) by 
the following Taylor expansion: 



where v is the barycentric frequency, v is the frequency derivative, and t is some ref- 
erence epoch. We were able to whiten our phase residuals with just a frequency deriva- 
tive (see Fig. El bottom). We find v = 442.33850(5) Hz, v = 3.7(2) x 1(T 6 Hz s" 1 at 
to=53535. 463215 MJD (UTC). The phase residuals after subtracting the model given by 
Eq. [2] are shown in Fig. [3] (bottom panel). Notice that the pulse train lasts for 576 s; this is 
quite long when compared to those of Type I X-ray bursts. We folded all the events within 
the 576 s pulse train using our best fit v and z>, and the resulting pulse profile is shown 
in Fig. HI Notice that the pulse profile is highly sinusoidal. Fitting a sinusoid to the pulse 
profile returned a reduced \ 2 value of 0.76 for 13 degrees of freedom, affirming the absence 
of any harmonic content. From the fit we obtain a fractional pulse amplitude of 2.1±0.1%. 

It is possible that the frequency drift we have observed is entirely due to orbital motion. 
The observed spin frequency of a neutron star in a circular orbit is given by 



where, u is the barycentric frequency at time t = t , i> ns sinz is the projected velocity of the 
neutron star, T is the orbital period, and <po is the orbital phase at time t = t Q . Unfortu- 
nately we do not have a long enough data set to place interesting constraints on T or v ns sin i; 
however, we can determine whether orbital modulation is consistent with the observed fre- 
quency drift. For example, if we assume, for solely demonstrative purposes, that this system 



<f>(t) = (f>{to) + v(to){t - t ) + -v(t )(t - t ) 2 + ■■■, 



(2) 




(3) 
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Fig. 3.— Pulse Frequency Evolution. Top: The burst lightcurve as shown in Fig. [TJ 
Middle: Dynamic Z 2 power spectrum. Notice that the pulsations remain highly significant 
(with Zf > 16) for 576 s. The solid line represents our best-fit frequency model consisting of 
just a frequency and a frequency derivative term. Bottom: Phase residuals after subtracting 
our best fit frequency model. Notice that the pulsations are highly coherent. 



- 8- 



36000 rr 



CO 
-t— ' 

cz 

o 
O 



35500 -ir 



35000 



34500 



34000 



ii if 



33500 L 1 1 1 1 

0.5 1 1.5 2 

Phase (Cycles) 

Fig. 4. — The pulse profile obtained by folding all events in the 576-s long pulse train (i.e., 
the interval where the pulsations were highly significant as determined from the Z 2 power 
spectrum) on the ephemeris determined in our timing analysis. The pulse shape is highly 
sinusoidal. The solid line represents the best-fit sinusoid. The fit yielded x 2= 0.76 for 13 
degrees of freedom, affirming the absence of any harmonic content. 
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has a projected velocity comparable to the one found by IStrohmayer fc Markwardtl (120021 ) 
for 4U 1636—53 (v ns sin i = 136 km _1 ), then fitting Eq. [3] to our frequency time series yields 
a good fit with a reasonable orbital period, i.e. T =7±1 hr. Hence, binary motion can 
account for the observed frequency drift. Observing another outburst would be essential in 
ultimately constraining the orbital parameters of this system. 



3.3. Spectral Analysis 



The timing analysis seems to suggest that this event is in fact very different from a 
canonical burst oscillation. To study it further we performed a spectral analysis. We analyzed 
the spectral evolution of the burst by breaking the burst up into ten intervals. The exposure 
of each interval was selected by demanding that each one contained an equal number of 
photons. Using the last interval as an estimator of the persistent emission, we fit each 
interval to a simple photoelectricall y absorbed b l ackbo dy model while holding the column 
density fixed to the value found by iPooley et al.l (120021 ) for the optical reddening. We find 
only subtle evidence for spectral softening, significant only at the 1-er level. 



Using the distance estimate to the cluster found by lOrtolani et al.l (119941 ). d = 8.5 ± 
0.4 kpc, we were able to calculate the luminosity of the burst. At the start of the observation 



the flux was ~2xl0 erg s , which is ^O.lLg^d, where Lem is the Eddington luminsosity 
for a neutron star. The burst lightcurve was well fit by an exponential with a decay timescale 
of 22 s. Unfortunately the observation only caugth the tail end of the burst-like event, we 
could not therefore characterize the peak luminosity of the burst. 



4. Discussion 

4.1. Super Burst? 

The timing properties of the pulsations are reminiscent of those observed during super- 
bursts. In particular the highly coherent pulsations and the long pulse trains. The superburst 
interpretation would be very exciting given the rarity of such events. However, we could not 
unambiguously determine this is the case given that we only caught the tail end of the burst. 
To study this phenomenon further we extrapolated our luminosity model determined in our 
spectral analysis before the start time of the observation. We find that the burst reaches 
the Eddington luminosity ~50 s before the start of the observation. Now, if we assume that 
we caught the tail of the outburst it is very surprising that it dropped from r^L Edd down 
to ~.01L#dd in ~77 s. This is much too fast for a superburst-like event. If, on the other 
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hand, this is a precursor events, such as those seen in superbursts, than it is unclear why 
the persistent flux did not rise back to ^L^dd shortly after as observed in other superbursts. 
Thus, the spectral analysis of this event suggests it is not a superburst. 



4.2. Intermittent Pulsar? 



The timing properties of the pulsations share all the hallmarks of those observed during 
superbursts, but the spectral analysis makes this interpretation unlikely. If this event is not 
a superburst, then could it be that these highly coherent pulsations h ave eminated from an 
accretion-powered millisecond pulsar? Recently iGalloway et al.l (120071 ) discovered "intermit- 
tent" pulsations from the accretion-powered millisecond pulsar (AMP) HETE J1900. 1—2455. 
They found that the properties of this pulsar differed from those of the other six known 
AMPs. For example, the pulsations were only present in the first f ew months of its outburst 
as opposed to other AMPs which show pulsations throughout (see IWijnands 112004 for a re- 
view). In addition, the pulse amplitude increased at three points in time that were almost co- 
incident with the times of thermonuclear bursts. The pulsed amplitude subsequently decaded 
after the bursts. Thus, it is plausible that we h ave observed a similar phenomenon. Following 



the th eoretical work of ICumming et al.l (120011 ) and lPayne &: Melatod (120061 ). IGalloway et al. 



(120071 ) suggest that the reason why the pulsations in HETE J1900. 1—2455 sometimes "switch 
off" is because the magnetic field is burried by accreted material, thus it cannot channel the 
flow which gives rise to the pulsations in the first place 

If the event we have discovered is another example of these intermittent pulsations from 
an unidentified source in NGC 6440, then a future outburst should reveal more long pulse 
trains, and if they happen to be correlated with the times of Type-I X-ray bursts then this 
would solidify the connection between this source and HETE J1900. 1—2455. 



5. Conclusions 

We discovered a 442 Hz pulsation during a burst from an X-ray source in NGC 6440. 
We could not establish the energetics of the burst as RXTE only caught the tail end of the 
event, but based on the timing properties of the pulsation we conclude that it was not a 
canonical burst oscillation. In particular we find that: The pulse train lasted for 576 s, as 
opposed to the ~10 s long pulse trains observed in Type I X-ray bursts. The signal was 
highly coherent and drifted down only by 2.1 x 10~ 3 Hz, as opposed to Type-I X-ray burst 
which exhibit drifts which are three orders of magnitude larger. The pulsations share all 
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the properties of those observed during superbursts; however, the complete energetics of this 
event make this interpretation implausible. We conclude that the long pulse train is most 
likely an "intermittent" pulsatio n from an accreting m illisecond pulsar such as those seen only 
in HETE J1900. 1-2455 thus far jGallowav et al.ll2007h . We not e that NGC 6440 is known to 
harbor a transient X-ray sourc e with a 409.7 Hz spi n frequency (IKaaret et al.ll2003l ) However, 
X-ray imaging of the cluster by lPooley et al.l (120021 ) has shown that NGC 6440 contains many 
X-ray sources. Thus, the burst presented here most likely emanated from a different object 
in the cluster. 
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